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Abstract. We discuss recent lattice results on in-medium properties of hadrons 
and focus on thermal properties of heavy quark bound states. We will clarify 
the relation between heavy quark free energies and potentials used to analyze 
the melting of heavy quark bound states. Furthermore, we present calculations 
of meson spectral functions which indicate that the charmonium ground states, 
J/ip and T] c , persist in the quark gluon plasma as well defined resonances with no 
significant change of their zero temperature masses at least up to T ~ 1.5T C . We 
also briefly comment on the current status of lattice calculations at non- vanishing 
baryon number density. 



1. Introduction 

Two quite different aspects of the studies of QCD thermodynamics on the lattice 
gained most attention in the heavy ion community during recent years - studies of 
the QCD phase diagram at non-zero baryon chemical potential and the analysis of 
the influence of a thermal heat bath of (quarks and) gluons on basic properties of 
hadrons, e.g. their masses and widths. In this review we will focus on the latter 
aspect of lattice studies but will also briefly discuss studies of the QCD phase diagram 
in this Introduction. 

1.1. Finite density QCD 

There has been considerable progress in extending studies of the QCD phase diagram 
to small but non-zero values of the baryon chemical potential /i^ QJQIISI- This led 
to a first analysis of the equation of state at non-vanishing baryon number density 
01 03 |S] and estimates for the location of the chiral critical point (second order phase 
transition) at which the transition to the high temperature phase of QCD turns from 
a rapid crossover to a first order phase transition. Recent calculations suggest a critical 
value [i% lt ~ 400 MeV IE] ■ This is much larger than the estimate for the baryon 
chemical potential which characterizes the chemical freeze-out of hadron resonances at 
RHIC, fi^ eeze ~ 29 MeV |H] and suggests that the transition at /is = as well as in 
the almost baryon free region produced at mid-rapidity at RHIC is characterized by a 
rapid crossover from a hadronic resonance gas to the quark gluon plasma. The analysis 
of the equation of state for fig > shows, on the one hand, that thermodynamics 
in the low temperature phase of QCD is well approximated by thermal properties 
of a hadronic resonance gas ^Hj- On the other hand, it also shows that the strong 
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increase of baryon number fluctuations, which for fig > and T < T c is consistent 
with the rise of baryon number fluctuations in a resonance gas, is suppressed once 
the system enters the plasma phase The resulting pronounced peak in the baryon 
number fluctuations indicates the proximity of the chiral critical point in the QCD 
phase diagram. Although it has been suggested since a long time that these kind 
of fluctuations should result also in visible event by event fluctuations in heavy ion 
collisions, it remains to be seen to what extent the experimentally observed increase 
in the K + /ir + ratio ^I] an d its fluctuation can be attributed to the QCD phase 
transition. 

1.2. In-medium properties of hadrons 

Confinement and chiral symmetry breaking are main characteristics of QCD which 
explain the basic properties of the hadron spectrum in the vacuum. Qualitative 
features like the existence of a light Goldstone particle, the pion, and the excitation 
spectrum of heavy quark bound states can be understood through the mere existence 
of a non- vanishing, spontaneously generated chiral condensate (tptp), and the linearly 
rising heavy quark potential with a non-vanishing string tension a, respectively. 
Quantitative properties of the spectrum depend, however, on the actual value and 
T-dependence of (ify) and a as well as on other properties of QCD like the running 
of the coupling with temperature or the generation of a thermal, Debye screening 
mass. It recently has been pointed out that despite of chiral symmetry restoration 
and deconfinement it is well conceivable that a complicated quasi-particle structure 
may persist in the QCD plasma phase at temperatures T > T c 13 . In fact, the 
non-perturbative structure of the plasma phase for T < 2T C is already evident from 
the strong deviation of bulk thermodynamic quantities , e.g. the pressure, from their 
perturbative form |14) . 

The influence of a thermal heat bath on hadron properties immediately becomes 
evident when one compares Euclidean time (r) correlation functions of mesons in 
different quantum number channels. At low temperature chiral symmetry breaking 
and the breaking of the axial U(l) symmetry lead to a splitting of hadronic states, 
which would be degenerate otherwise and would then also have identical correlation 
functions. In Figure ^ we show pseudo-scalar (tt) and scalar (5) meson correlation 
functions in a gluonic heat bath (quenched QCD) at temperatures below (left) and 
above (right) the deconfinement transition temperature; 



with Jh denoting a hadron current with the appropriate quantum numbers of the 
hadronic channel H . The spectrum in the pseudo-scalar and scalar channels differs at 
low temperature because of the explicit breaking of the axial £7(1) symmetry which 
consequently leads to quite different correlation functions in Figure^left). The almost 
perfect degeneracy of both correlation functions above T c , however, suggests that this 
symmetry is effectively restored in the deconfined phase of QCD (Figure fright)). 

The hadron correlation functions, G#(t, T), are directly related to spectral 
functions <Jh(w, T) which contain all the information on thermal modifications of the 
hadron spectrum in certain channels of hadron quantum numbers, H, 
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G h (t,T)= [ 
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Figure 1. Scalar and pseudo-scalar correlation functions at T = 0.6T C (left) 
and T = 1.5T C (right). Shown are results from calculations with light quarks on 
a 64 3 X 16 lattice in quenched QCD. 



Lattice studies of in- medium properties of hadrons greatly advanced in their predictive 
power through the exploitation of the Maximum Entropy Method (MEM) [T^. This 
allows the reconstruction of o~h{u,T) at non-zero temperature for light as well as 
heavy quark bound states. We will concentrate here on the heavy quark sector and, 
in particular, will discuss in how far potential model calculations, which generally 
led to an early dissociation of heavy quark bound states ]TQ can be made consistent 
with recent studies of charmonium spectral functions in the high temperature phase 
of QCD |17l 1181 [T5] . which find charmonium bound states even at T ~ 1.5 T c . This 
will require to reconsider the relation between heavy quark free energies calculated on 
the lattice and temperature dependent potentials used in potential model calculations. 
We will start with a discussion of this problem in the next section and present recent 
results on charmonium spectral functions in section 3. 



2. Color averaged heavy quark free energies and color singlet potentials 

The original idea of charmonium suppression as a signature for plasma formation in 
heavy ion collisions 20 relies on the occurrence of fundamental changes in the heavy 
quark potential (vanishing string tension, Debye screening) in the high temperature 
phase of QCD. While there is no doubt that this does happen in QCD it was clear 
right from the beginning of studies of ■//'(/'-suppression in potential models |2J that 
a quantitative determination of dissociation temperatures depends on the detailed 
structure of the heavy quark potential at high temperature. Information on the latter 
can be extracted from an analysis of heavy quark free energies |22| calculated on 
the lattice which describe the change in free energy of a thermal medium due to the 
presence of static quark and anti-quark sources, 

e - P ^T ) = ^V)^l_ ^ (3) 

x,y;\x— y\— r 

Here the static sources are represented by Polyakov loops L{x) and L'(y), and 
Zq q (T,V), Z(T,V) denote the QCD partition functions in presence or absence of 
these sources, respectively Quite often Fq q {r,T) itself has been interpreted as the 
heavy quark potential at finite temperature and as such has been used in potential 




Figure 2. The color averaged heavy quark free energy at several values of 
the temperature (left) and the temperature dependence of its asymptotic value 
at large distances (right): The left hand figure shows results for T/T c = 
0.9, 0.94, 0.98, 1.05, 1.2, 1.5 (from top to bottom) obtained in quenched QCD. 



model calculations. This, however, ignores the fact that aside from modifying the 
energy of the thermal heat bath the external sources also change the entropy of the 
system. The change of energy and entropy will depend on temperature as well as the 
separation of the gg-pair, 

The relative importance of energy and entropy contributions to the free energy at 
different qq separations becomes evident after a proper renormalization of Polyakov 
loops |23 EH When renormalizing L(x) by matching F qq (r,T) = F qq (r, T) - T In 9 at 
short distances to the zero temperature heavy quark potential [231 one readily finds 
that at short distances the temperature dependence of F is weak. This indicates that 
after subtraction of a simple statistical entropy term (T In 9) the short distance part 
of the free energy reflects the change in energy due to the presence of static quark anti- 
quark sources. This, however, is not sufficient to obtain the (potential) energy at all 
distances. At large distances a strong temperature dependence remains even after the 
subtraction of Tin 9. In fact, F qq (r,T) and F qq (r,T) will approach (— oo) for T — > oo 
(see Figure ■ As seen in FigureEb, f 00 (T) = F qq (oo,T) decreases almost linearly 
in T, i.e. ^^(T) ~ — 0.2T. This shows that the free energies contain a substantial 
entropy contribution at large distances while the energy contribution almost vanishes, 
which according to Equation (0J is reflected by the weak temperature dependence of 
Foo(T)/T. This complicated r-dependent entropy contribution makes a direct use of 
free energies in potential models questionable. 

A further obstacle for a direct use of lattice calculations of heavy quark free 
energies in potential models arises because at non zero temperature F qq does not 
directly give the singlet free energy from which one could determine the singlet energy 
or potential used in phenomenological approaches. The free energy defined through 
the Polyakov loop correlation function in Equation Q gives a weighted average over 
singlet and octet contributions [22] ■ At non-zero temperature the definition of a singlet 
free energy requires gauge fixing. At present most studies of the singlet free energy 



News from Lattice QCD 



5 




0.5 1 1.5 2 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



Figure 3. Singlet free energy and heavy quark energy in quenched QCD 
calculated in Coulomb gauge at T = 1.5 T c (left) and heavy quark energy at 
several other values of the temperature (right). 

have been performed in Coulomb gauge ||, where the singlet free energy is given by, 
e-^> = i £ l(TrL(x)L\y)) . (5) 

\x— y | — r 

Some results for T = 1.5 T c are shown in Figure |31 As expected from the discussion 
given above it is clear that at short distances free energy and energy will coincide while 
at large distances the asymptotic value of the heavy quark energy will be larger than 
the corresponding free energy value as the entropy contribution has been removed. 
This suggests that heavy quark bound states could survive in the QCD plasma up to 
temperatures significantly higher than those estimated previously in potential model 
calculations, which have been based on input obtained from calculations of the heavy 
quark free energy. To quantify these statements it is necessary to extent the analysis 
discussed here in the context of quenched QCD also to the case of QCD with a realistic 
light quark spectrum. 

3. Spectral analysis of hadronic correlators at high temperature 

It should have become clear from the discussion in the previous section that an analysis 
of thermal properties of heavy quark bound states in terms of potential models and 
probably even the determination of singlet potentials from lattice calculations will 
require additional phenomenological consideration. On the other hand, as pointed 
out in the Introduction, an ab-initio approach to thermal hadron properties through 
lattice calculations exists and is based on the calculation of thermal hadron correlation 
functions, G#(t, T). As indicated in Equation @ Gh(t,T) is related to the spectral 
functions ajj(uj,T). However, an inversion of this integral equation is generally 
not possible because lattice calculations only yield information on Gh (t, T) at a 
finite, discrete set of Euclidean time steps, t^T = k/N T with k = 0,1,... N T — 1, 
with N T denoting the temporal extent of the lattice. It is, however, possible to 

| To avoid an explicitly gauge dependent formulation one could introduce additional spatial Wilson 
lines to connect the quark anti-quark sources (cyclic Wilson loop). This trades "gauge dependence" 
for "path dependence" 1251 . However, also the singlet free energy calculated in Coulomb gauge can 
be understood in terms of an appropriately defined gauge invariant operator 1261 . 
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Figure 4. Thermal correlation functions above T c in axial vector (AX) scalar 
(SC), pseudo-scalar (PS) and vector (VC) channels normalized to reconstructed 
correlation functions which are based on spectral functions calculated at T* = 
0.75T C (left) and T* = 0.9T C (right). Note the different scales on the ordinates of 
both figures. 



determine the most probable spectral function which describes the calculated data set 
{Gir(Tfc, T) | k = 0, N T — 1} and respects known constraints on o~h (w, T) (positivity, 
asymptotic behavior, ...). This can be achieved using a Bayesian data analysis, e.g. 
the Maximum Entropy Method (MEM) 27 . In the context of QCD calculations this 
has been introduced in |15|. 

A first feeling for the influence of a thermal heat bath on the structure of 
correlation functions and the size of medium modifications of a thermal spectral 
function can be obtained by comparing directly the numerically calculated correlation 
functions at temperature T with a correlation function, 67 rccon (rfc, T), constructed 
from the spectral function calculated at a smaller (zero) temperature, o-h{u,T*) with 
T* < T, i.e. 

G recon {T,T)= duja„(uj,T) . r , . . ■ 6 

J smh(w/2T) 

In this way the trivial temperature dependence of the integration kernel is taken care 
of and any remaining discrepancies between the reconstructed correlation function 
and the actually calculated correlation function at temperature T can be attributed 
to changes in the spectral function. Some results for charmonium correlation functions 
are shown in Figure 01 • 

The ratio of correlation functions shown in Figure 0] would equal unity, if the 
spectral functions would not depend on temperature. However, as can be seen the 
correlation functions for P-state charmonia (left) show strong thermal modifications 
already at temperatures slightly above T c . On the other hand, the 5-state correlators 
(right) show only little modifications up to T ~ 1.5 T c . In fact, no temperature 
dependence is visible in the pseudo-scalar channel (rj c ) up to T = 1.5 T c while 
modifications in the vector channel ( J/4>) are of the order of 10% at this temperature. 

The reconstruction of charmonium spectral functions has been performed using 
technically different approaches, e.g. using point-like ^| El or smeared (2H1 hadron 
sources on isotropic ^H] or anisotropic |171 IT^l lattice with standard Wilson fermion 
(17) or improved Wilson fermion [181119} actions. The different calculations agree to the 
extent that no significant modification of 5-state spectral functions is observed up to 
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Figure 5. Thermal vector spectral functions (J/tp) from MEM analyses of meson 
correlation functions calculated in quenched QCD on anisotropic |17| (left) and 
isotropic 1181 (right) lattices. The difference in scale on the ordinates of both 
figures is due to different normalizations used for the hadronic currents. 



1.5 T c , i.e. the J/tp survives as a narrow bound state with unchanged mass up to this 
temperature. This is shown in the upper part of Figure[SJ Current lattice calculations^ 
differ, however, on the structure of spectral functions for larger temperatures (lower 
part of Figure |3J) ■ While it is concluded in [T7| that the J/tp resonance disappears 
quite abruptly at T ~ 1.9 T c the analysis of ^B] suggests that the resonance disappears 
gradually; a resonance peak with reduced strength is still visible at T = 2.25 T c and 
finally disappears completely at T — 3 T c . 

To get control over the detailed pattern of dissolution of the heavy quark 
resonances clearly requires more refined studies. It also should be noted that so far all 
existing lattice studies have been performed in the quenched approximation. Although 
virtual quark loops are not expected to modify the qualitative picture obtained from 
these calculations in the heavy quark sector of QCD one clearly has to understand 
their influence on a quantitative level. 

4. Conclusions 

Studies of spectral functions of charmonium states suggest that the cc ground states, 
J/tp and rj c , exist as well localized resonances with masses essentially equal to their 

% Note that current lattice calculations of spectral functions are strongly influenced by lattice cut-off 
effects which show up most strongly at large energies, u. In fact, only the first, low energy peak in 
the spectral functions shown in Figure [S] is physical and insensitive to changes of the lattice cut-off. 
The other two peaks have been shown to be lattice artifacts arising from "Wilson doublers" I2i)l . 
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vacuum values at least up to temperatures T ~ 1.5 T c . Whether these states 
disappear abruptly at higher temperatures or whether their contribution to the 
spectral functions gradually weakens requires further studies. 

We also discussed the reanalysis of heavy quark free energies and showed that 
these cannot directly be identified with the potential energy of a static quark anti- 
quark pair. The elimination of entropy contributions to the free energy leads to 
(potential) energies which generally have a larger dissociation energy for a cc-pair. 
Taking this into account in potential model calculations, which in the past suggested 
a rather early dissociation of heavy quark bound states, can well lead to dissociation 
temperatures consistent with results obtained from the analysis of spectral functions. 
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